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Abstract 

The renormahzation of the effective QCD-Hamihonian theory for the quark- 
antiquark channel is performed in terms of a renormahzed or fixed-point 
Hamiltonian that leads to subtracted dynamical equations. The fixed point- 
Hamiltonian brings the renormalization conditions as well as the counterterms 
that render the theory finite. The approach is renormalization group invari- 
ant. The parameters of the renormalized effective QCD-Hamiltonian comes 
from the pion mass and radius, for a given constituent quark mass. The Is 
and excited 2s states of uq are calculated as a function of the mass of the 
quark q being s, c or b, and compared to the experimental values. 
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I. INTRODUCTION 



The effective mass operator equation for the lowest Light-Front Fock-state component 
of a bound system of a constituent quark and antiquark of masses mi and m2, obtained in 
the effective one-gluon-exchange interaction approximation, which is our starting point, was 
reviewed in Ref. [|T|. The breakthrough to simphfy this equation in the spin zero channel 
was achieved in [§] by formulating the t|-model, which is reduced to 



X 1 — X 



ij{x,k_i 



f dx'dk\9(x')9(l — x') /4mim2 a , \ , / / r/ \ 

where M is the mass of the bound-state and is the projection of the light-front wave- 
function in the quark- antiquark Fock-state. The mean four-momentum transfer is 0] 

Q\x,k^,x',k'^) = -^[ik^-k[f-{k2-k',f] . (2) 

The coupling constant a for the Coulomb-like potential and A is the bare coupling constant of 
the Dirac-delta hyperfine interaction. For the purposes of this work, we leave out the energy 
transfer in Eq.(|^), that leads to aditional singularities for the tridimensional momentum 
going to infinity, that we are not going to treat here. The effective mass operator of Eq. (|l]) 
acting in the quark-antiquark sector, has also been obtained by iterated resolvents, which 
were derived in |Q and presented in greater detail in 0, that allows to express systematically 
the higher Fock-state components of the wave-function in functionals of the lower ones. In 
this way, the higher Fock-state components can be retrieved from the qq projection, and 
the full complexity of the QCD theory is in principle described by the effective Hamiltonian 
acting in the lowest Fock-state component [Q. 

For convenience, the Sawicki tranformation, first derived for equal masses and consis- 
tently formulated for unequal masses in Ref. [Q, is applied to Eq.(|I|) which allows to rewrite 
it in the instant-form momentum basis. It is useful in this case since the momentum transfer 
is approximated by a rotational invariant form given by Eq.(0): 



.(fe) = , (3) 



and the Jacobian of the transformation of (x, /c^) to A; is: 



dxdk^ = ^^^-rj^dk , (4) 
mrA{k) 



with the dimensionless function 



Mk) = —1^^ , (5) 



and the reduced mass rrir = mim2/ {nii + 7712). The individual energies are Ei = ymf + k"^ 
(z=l,2) and k = \k\. 

The mass operator equation in instant form momentum variables is given by: 

MMk) = [E, + E^f m - I dk' I ^ , " , ^ I Vik') , (6) 



where = mi + m2, the phase-space factor is included in the factor 1/ ^ A{k)A{k') and 
^x{1-x)^Ij{x, k±) = ^A{k)ip{k). 

The mass operator equation needs to be regularized and renormalized in order to 
give physical results. In Ref. the delta function was smeared out to a Yukawa form, 
regularizing Eq.(P) and the parameters found from the pion mass and radius. In principle 
the size parameter in momentum space of the Yukawa potential should be let to infinity while 
the physical input is kept constant. To make vanish the dependence on the size parameter 
is a non trivial task, which is the purpose of this work. Here we obtain a renormalized 
form of the equation for the bound state mass, which is i) invariant under renormalization 
group transformations, ii) the physical input is given by the pion mass and radius in) no 
regularization parameter. 

We are going to apply the recent renormalization techniques developed in the context 
non-relativistic Hamiltonian theory to Eq.(^, since it defines an effective Hamiltonian 
for the quark- ant iquark dynamics. The renormalized T-matrix is the solution of a subtracted 
scattering equation, which the physical input is given by the T-matrix at some reference scale 



/i. The scheme is invariant under the change of the arbitray scale /i and consequently the 
inhomogeneous term of the subtracted scattering equation satisfies a renormalization group 
equation, which expresses the matching of the theories at the scales /i and fi + dfi 0. We 
will find the mass of the bound state from the pole of the renormalized scattering matrix 
defined from the mass operator Eq.(^. 

This work is organized as follows. In sec. II, the operators for the Coulomb-like and for the 
singular interactions appearing in the mass operator are defined, and the Lippman-Schwinger 
equation for the T-matrix related to the given mass operator is written. In sec. Ill, we show 
how to renormalized the effective theory defined by the mass operator in sec. II through the 
definition of a renormalized or fixed point singular interaction. This procedure is equivalent 
to use subtracted scattering equations. We discuss the renormalization group invariance of 
the method. The explicit form of the renormalized T-matriz is obtained (appendix A) and 
the physical input is the pion mass which determines the excited states as well as the mass 
of the other mesons. The physical observables of the renormalized effective theory does not 
depend on the subtraction point. In sect. IV, the calculation of the pion charge radius is 
discussed and although we have simplified the spin dependence in the dynamical equation 
it is important in the evaluation of the radius. For this purpose we have used an effective 
pseudo-scalar Lagrangian to construct the spin part of the pion wave-function, and turning 
off the Coulomb-like interaction we retrieve a well knwon result valid in the soft pion limit 



14 1 . To gain insight we also write down the expressions where the quark spin is neglected. 



The numerical results are presented in sec. V. We have solved a nonrelativistic example for 
a Coulomb plus a Dirac-delta interaction compared to a model in which the Dirac-delta is 
substituted by a Yukawa. We also show how this effective theory, including a Dirac-delta, 
can mimic a finite range theory, calibrated to the pion mass. Then, we present results for 
the pion charge radius and for the mass of the excited state, that we idendified with the 
isovector vector mesons. Finally, in sec. VI, we present our conclusions. 
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II. DEFINITIONS 



For our puposes it is convenient to work in an operator form of Eq.(^): 

[mI + V + V^)\^>=M^\^> , (7) 

and the free mass operator Mq is the sum of the energies of quark 1 and 2, V the Coulomb- 
hke potential and is the Dirac-delta interaction in the non-relativistic hmit. The matrix 
elements of these operators are given by: 

< k\V\k' > = -^""j , " , , (8) 

and the short-range singular interaction 

<k\V^\k'>=<k\x> — <x\k'>=— — — , ^ =■ (9) 

rrir ^A{k) ^A{k') 

The phase-space factor A{k) is defined by Eq.(^, the square momentum transfer Q"^ comes 
from Eq.(l^). For convenience in the formal manipulations of the next section, the form-factor 



of the separable singular interaction is introduced and defined by < k\x >= l/yA(k). 
The T-matrix is obtained from the Lippman-Schwinger equation: 

T(M2) = V + + {V + V^)Gi^\M'^)T{M^) , (10) 

where M is the mass of the scattering state and the Green's function with outgoing wave 
boundary condition is 

In the next section we will obtain a renormalized form of T{M'^). 

III. RENORMALIZATION OF THE EFFECTIVE THEORY 



The ideas that have been developed in Refs. 0^ to construct a renormalized T-matrix 
in non-relativistic Hamiltonian theory, can be applied in the case where the two-body inter- 
action can be splited in a regular potential of finite range iV) and a Dirac-delta singularity 



as in Eg. ([T0|) . In this case, the physical information at the subtraction point will be intro- 
duced through the renormalized mass operator as well as all the counterterms that render 
finite the T-matrix equation (plQl). The physical information at the subtraction point is the 
T-matrix corresponding to the Dirac-delta interaction. This renormalization approach has 
been applied in the two-nucleon system to calculate the T-matrix of the one-pion-exchange 
potential, without the necessity of regularization or form-factors 0. 
The renormalized interaction is given by: 

Vn = V + Vl, (12) 

where the renormalized Dirac-delta interaction is given by: 

1 

l + r^(/J)G(+)(/i2) 



(/i 
1 



n=0 



1 + G[,+H/.2)T^(/.2) 

oo 

T^{y^')Y.[-Go{^')T^{^')Y ' (13) 



and T^(/i ) is the renormalized T-matrix of the Dirac-delta interaction, with matrix elements 
given by: 

<p\Ti{l/)\q>=<p\X > A7^(/i') < X\q> , (14) 

where X-jil^"^) is the renormalized strength of the Dirac-delta interaction at the mass scale 
/i^. In the non-relativistic limit the form-factor x(g) = 1, and the renormalized interaction 
becomes the Dirac-delta. The scattering equation with the renormalized interaction appears 
in a subtracted form [^-0 after a little rearrangement of terms, in which all the divergent 
momentum integrals are removed and it is written as 

TniM') = Tnifi') + Tn{^') {G^^\m') - gJ+)(/x^)) Tn{M^). (15) 

It is presented here in a general way, where we have dropped the 5 superscript, just to remind 
the reader that for a regular potential Eq.(P^ is completly equivalent to the traditional 
Lippman-Schwinger scattering equation. 



The renormalized interaction is independent of the subtraction point, i.e., the physics 
expressed by the renormahzed interaction is invariant by changes in the arbitrary renormal- 
ization point, this physical requirement is given by: 

d 



^V-O; (16) 

quahfying the interaction as the fixed-point of Eg . (|16|) , which imphes that the T-matrix 
found from the solution of 

Tn{M^) = V + Vl+(y + G'C\M^)Tn{M^)) , (17) 

is invariant under dislocations of the subtraction point. Consequently, the renormalized 
coupling constant of the Dirac-delta interaction changes as the subtraction point moves, 
according to the Callan-Symanzik equation 

^Tiin') = -Ti{^^')-—-^-—-,T^{^^') , (18) 



obtained from Eq.(^). 

Although the renormalized interaction is not well defined for singular interactions, the 
resulting T-matrix obtained by solving Eq.(|T^) is finite. This gives a posteriori justification 
for the formal manipulations used in (|13D and (plTf). The sum in the expression of the renor- 
malized interaction (|13D explicits all the counterterms which exactly cancels the infinities 
in the momentum integrals of the scattering equation (0), while introducting the physical 
information through the value of the renormalized strength of the Dirac-delta interaction. 
We choose A^(/x^), in the following, in accordance with the physical value of the pion mass. 
We observe that, instead of working formally with the operator V^, we could use an ul- 
traviolet momentum cut-off (A), by defining in this way a regularized interaction. After 
the construction of the T-matrix regularized equation one could perform the limit A oo, 
ariving at the same results as the ones obtained directly with the use of the renormalized 
interaction. 
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The solution of the scattering equation (|T^ is found by using the two potential formula 
in terms of the T-matrix of the regular potential V, T^(M^) and the renormalized T-matrix 
of the Dirac-Delta interaction, see Appendix A, which results in 

{l + T^{M^)Gi^\M^)) \x >< x\ {Gi^\M^)T^{M^) + l) 

A^^(/.2)- < x\ (iPTihu- - ^ITOw-) \x>-< x|gW(M2)T^(M2)gW(M2)|x > ' 

(19) 

where the T^(M^) is the solution of the Lippman-Schwinger equation 

T^(M2) = V + VGi^\M^)T^{M^) . (20) 

The structure of Eq. ([T9|) allows one more subtraction in the denominator, which turns 
faster the convergence of the momentum integral in the term where T^(M^) is present. This 
subtraction is appropriate if the potential V has a Coulomb or Yukawa form, 

<p\v\q>= 2^,1 (21) 

+ |p — q] 

Thus, we define the renormalized strength of the Dirac-delta interaction at the subtraction 
point such that, 

A^'(^') = \n\fi')+ < x\Gi-'\fi')T''ifi')Gl;\fi')\x > ; (22) 



and introduce the physical information in the renormalized T-matrix (p!9|), through the value 
of 

Substituting Eq . (|2^) in (p^Of) , we obtain the renormalized T-matrix written as: 



+ (i + T''{M^)G^^\m^)] \x > tn{M^) < x\ (G^^\m^)T'' {M^) + l) , (23) 



where, 

tn\M') = Xnif^')- < X\ {G^'^^Hm') - G^W(/i^)) \x > , (24) 



and the interacting Green's function is 

G^(+)(M2) = G^^\m^) + G^^\m^)T'^ {M^)G^^\m^) . (25) 

We use the renormahzation condition that at the pion mass, M = rriT^, the T- matrix, for 
mi = 1712 and mi = niu = rrid has the bound-state pole, consequently 

t^^(m^)=0, (26) 

and choosing /i = m^^ for convenience, which implies that 

'^n\0 = 0- (27) 



The invariance of the renormalized T-matrix ( p3D under the dislocation of the subtraction 
point, just reads that 

d 



and from 



MM') = , (28) 



Tnif^") = < X\ (G^W(/.'^) - G^W(^^)) \x > . (29) 



At the general subtraction point t-j^^fi^) = X-j^^ji^) and the renormalized T-matrix at jj? 
is given by: 

(l + T^(/i2)G{,+)(^^)) \x > V(/i') < x\ + l) • (30) 

The full renormalized interaction can be written in a form analogous to Eg . (|T^) : 

oo 

Vn = Tnif^') E [-G^o'\f^')Tn{f^')Y , (31) 

n=0 

and with that, one could obtain the equation for the renormalized T-matrix (|1^) in the 
subtracted form (|TB|), displayed again below 
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We observe that one could equally well construct the Callan-Symanzik equation for 



-^TTi{fx^), by performing the limit of M — /x in ([T5|), finding 



d 



(32) 

gives the dependence 



with the boundary condition given by Eq. (pOD . The solution of Eq. 
of At^ on the subtraction point /i as expressed by Eq. p9|) . 

Now comes a subtle point: It is important to realize that the renormalization condi- 
tion given by A^^(m^) = 0, and considering Eq. p^) , the bare strength of the Dirac-delta 
interaction is given by 

mrX,L = [< xlG^'^^HrnDlx >], ^ (33) 

where the Green's function of the Coulomb like interaction is calculated for equal mass 
constituent quark and antiquark, i.e., nii = 1712 = rriu = m^. The bare coupling constant 
in the form explicited by (^) is sufficient to render finite the T-matrices from Eq.(^). 
calculated with different quark masses. 

The bound-state mass, Mf,, of quark-antiquark states with also different quark masses 
are the position of the poles of the renormalized T-matrix Eq. (p3D , which implies 



rrir 



<x\G''{ml)\x> 



{mu,m-) 



rrir 



<x\G''{Mi)\x> 







(mi ,7712) 



(34) 



The zeroes of give the position of the zero angular momentum states of the bound quark- 
antiquark systems, for different quark masses and excitation. It is easy to imagine that in 
the vicinity of a bound-state of the Coulomb-like potential, because of the presence of the 
pole in , the function is rapidly varying and it is infinity at the Coulomb bound state 
and because of the change in sign of it will necessarily presents a zero, and a bound 
state of the whole potential. 



IV. PION CHARGE RADIUS 



The renormalized effective theory defined by Eqs.( PB| , |2^ ) with the bound state mass 
equation (|3^) , once the pion mass is known has only one free interaction parameter, a, the 
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strength of the effective one-gluon-exchange potential. The value of a can be found from the 
pion charge radius. The pion wave-function in the effective theory comes from the residue 
of the T-matrix, Eg. (pSf) , at the pion pole, such that 



where 



^^^^ Mk), (35) 



^ x(l — X 



=< k\Go{ml) (l + T^'imDGoiml)) \x > , (36) 

and 

4k,) = 1 + ^^ (37) 

from Eq.(|^). The absolute normalization of the qq Fock-component of the pion wave-function 
(^) is such that the asymptotic form is given by the first term of (^), with residue at the 
pion pole equal one. We impose such normalization condition to be consistent with the soft 
pion limit (m^=0) for the electromagnetic form factor when the Coulomb type interaction 
goes to zero, as it will be shown below. 



A. Including Quark Spin 

The pion electromagnetic form-factor is obtained from the impulse approximation of the 
plus component of the current (j^ = j° + j^) in the Breit-frame with momentum transfer 
g"*" = and = —(f. The leptonic decay constant of 7r+ (/^) is a physical quantity 
which depends directly on the probability to find the quark-antiquark Fock-state in the 
pion wave-function and consequently properly normalizes it once the empirical value of 
is given. Computing the pion form-factor from an effective Lagrangian, described below, 
the value of gives the normalization of the form-factor. In this case the qq component of 
the pion wave-function is normalized such that, in the vanishing limit of the Coulomb- type 
interaction, it retrieves the asymptotic form Eq.(|^) |12|. The coupling of the pion field to 
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the quark field are taken from an effective Lagrangian witli pseudo-scalar coupling with the 
pion quark coupling constant given by the Goldberger-Treiman |11[] relation at the quark 
level, 

^eff = -iY^.qj^rq , (38) 

In 

and /,r= 93 MeV is the pion weak decay constant. 

The general structure of the qq bound state forming the pion comes from the pseudo- 
scalar coupling, and we will use such spin structure in the computation of the Feynman 
triangle diagram which expresses the impulse approximation to compute the pion electro- 
magnetic current. 



7^ 



{k — p'^y — m? + ie {k — p^)^ — m? + ie ' 
where F^^lq'^) is the pion electromagnetic form-factor, the pion momentum in the initial and 
final states are defined by p° = p'^ and p'^j_ = —pn± = ^ in the Breit-Frame. iVc = 3 is the 
number of colors. 

The -|- component of the current is calculated from (pQl). It is chosen because after 
the integration over k- = k° - k^ the supression of the pair diagram is maximal for this 
component in the frame where g+ = and just the valence wave-function enters in the 
form-factor [|1^,|1^]. Although we are going to compute the integration in the - component 



of the moment assuming a constant vertex, one can identify in the expression how the wave- 
function correspondent to the non-constant vertex of Eq. (|35D should be introduced in the 



expression. The result is 

where the momentum fraction x = {p^ — k^)/p'^. The relative transverse momentum is 
given by: 

K± = {1- x){p± - k±) - xk_L, (41) 
12 



and K'^ = K± + (1 — The free mass operator for the q — q is written in terms of the 

momentum fraction and the relative perpendicular momentum 

MS = (42) 

and Mq is written as a function of K'j_. The expression for the pion form-factor gives 
the standard Drell-Yan formula once the bound-state wave-function of the constant vertex 
model, the asymptotic form, is recognized 

^^{x, K^) = . ^ , (43) 



which is the first term of the pion wave-function in Eq. (pG]). The second term in (|36D comes 
from the Coulomb-like potential, and implies in a contribution from the effective theory 
to the pion radius. The other factors in Eq.(^OD comes from the Melosh rotations of the 
individual spin wave-function of the quarks. 

The pion charge radius from Eq.(^) in the soft-pion limit with constant vertex, cor- 
responding to the wave-function of Eq.(|^) with = 0, gives the well known result of 

1 

= [6^F^(g2)|g2^o]^ = v^/(27r/^) from Ref. H. The form factor (ID in the soft-pion 



limit with constant vertex, for q = 0, reduces to expression for |T2|, 

U = TTT^.^N, / , / d'KMx,K^) , (44) 



obtained from the computation of the leptonic decay transition amplitude of 7r+ with the 
effective Lagrangian (pSf ) and the wave function given by Eq. (|43|) . 

B. Neglecting Quark Spin 

In order to gain insight on the importance of the inclusion of the quark spin in the 
computation of the form-factor, we have simplified the numerator of Eq . (^0|) , which is the 
result from the Dirac algebra, by taking the limit of m ^ oo. In this way the quark 
spin is neglected and the form-factor reduces to the formula found for scalar particles, and 
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only the overall normalization, which depends on /^r is mantained. The wave-function is in 
the form expressed by Eq.(^) with A[k) = 1. With the above approximations the pion 
eletromagnetic form-factor is 

1 f 

FM^) = ^^N, / dx / d''K^Mx,K^)^P^{x,K'^) , (45) 

TT-^ Jo J 

with the relative transverse momentum 

= (1 - x){pi_ - k±) - xfcx, (46) 

and K'j^ = K± + (1 — x)q±. In case of the qq wave-function being supossed to be the complete 
pion wave-function, the form factor is written as |l|] 

F^{q^)=Af£dx J d^K^Mx,K^)Mx,K'_^) , (47) 

and from -^^(0) = 1 the normalization A/" is determined. 

The difference between the form factors defined by Eq. (^) and Eq. ( ^71) are the normaliza- 
tions of the qq Fock- component of the pion wave-function. In the first case the normalization 
is defined by while for the second case the qq Fock-state component is normalized such 
that Ft^{0) = 1. As we will see, Eq.(^) gives a too small pion radius, while from Eq.(^5D 
or Eq.(PD|) the pion radius can be described reasonably. The absolute normalization of the 
form-factor computed with the inclusion of the quark spin as expressed by Eq.(^0]) is unde- 
termined, since the integral diverges. However, the pion radius is finite with the magnitude 
of this Fock- component of the pion wave-function known from the empirical value of /^r, 
number of colors and the attributted constituent quark mass. 

To close this section, we observe that our aim is to fit the strength a of the Coulomb- 
like potential using the pion charge radius. For this purpose, we introduce the pion wave- 
function from Eq. ( pS]) in the form-factor expression, Eq. ( ^0]) and calculate the charge radius. 
Although the form-factor diverges, the charge radius is finite and for = and a = it 
retrieves the soft-pion limit. The pion mass and the effective Coulomb- like interaction give a 
correction to the soft pion limit, r^°-^*=0.58 fm towards the experimental result of 0.67±0.02 
fm . 
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V. NUMERICAL RESULTS 



A. Test Case 

We begin the section on the numerical calculations comparing our results for the s-wave 
bound-state energies which are consistent with []TU[ for the non-relativistic Coulomb plus 
repulsive Yukawa model 



[e - = - / dp'- 

n Jo p 



In 



(P - P') 



i\2 



In 



{p -p'Y + r] 



(48) 



{p + p'Y {p + p'y + 'rf_ 

In our numerical procedure we checked the results with up to 200 Gaussian-Legendre quadra- 
ture points and the interval —l<z<l was transformed to < A; < oo through the 
variable transformation k = c(l — z)/{l + z) with c about 1 . 

In Figure 1, we show our results for the first excited state e^^-* as a function of the ground 
state e^^-* and compare with the calculation of |10[ for 77 = 0.1 , 1 and 10. For these values 
of 1] we present results in Table I. Our precision is about 0.5%, that will suffice for our 
purposes. Also the calculation with the renormalized Coulomb plus Dirac-delta is shown in 
Figure 1 and Table I. The calculation with non-relativistic renormalized model corresponds 
to find the zeroes of the non-relativistic form of Eq. (|3^) , 

(49) 



J dqdp< q\ (G'^"'-(-e(")) - G^"'-(-e(i))) \p>= , 



with the non-relativistic resolvent G^^^{e) for negative energies obtained from the solution 
of 



G^™'(e) = G'o'(e) + G'^''(e)VbG^"^'(e) 

where < plVclg >= [n'^{\p — g|^)] ^ and the free resolvent G'o'^'(e) = [e — k 
In momentum space Eq.(^9|) is given by 

1 1 



21-1 



An / dp p^ 
Jo 



-An) 



.(1) _ p2 



POO POO 

+ Svr^ / dq q"^ dp p^ 
Jo Jo 



tc{p,q;-e(^^) 



tc{p,q]-e^^^) 



(_e(n) _ p2) (_g(n) _ g2) (.^(l) _ ^2) _ ^2) 



(50) 



0, (51) 
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and s-wave projected T-matrix in Eq.(^) is the solution of 

tc{p,q]e) = ——\n———^ + - dp -In —. 52) 

ZTr"^ pq [p + q) tt jo p [p + p ) e — p''^ 

We solve Eq. for n = 2 for each given e^^-* . Now, it is clear the motivation for plotting 

Figure 1 the binding energy of the excited state against the ground state energy. In Table 

I as well as in Figure 1, the renormalized model reproduce with less then 10% accuracy the 

model results of Eq . (^81) . In the limit of ?7 — oo, the renormalized model should work better, 

however in the test case of Eq. (|48|) this limit corresponds to a vanishing Yukawa interaction, 

and the results tends trivially the Coulomb value of the first excited state of 0.25 for the 

ground state value of 1, what in fact is observed in Figure 1 and Table I. 



B. Effective Pion Model 

The effective model of corresponds to use the non-relativistic phase-space A{k) = 1 
in Eq. (^) and a smeared delta-interaction of a Yukawa form: 



m'Mk) = Urn' + Ak'] ^ - A« / — ( + . !l r..J ^i^') , (53) 



Stt^ J m \{^k-k'Y r]^ + {k-k'y 
which was solved with parameters adjusted to fit the pion mass and the rho- meson mass, 
resulting m = 406 MeV, a = 0.6904, and r] = 1330 MeV f^. The value of t] was found from 
the condition that the first excited bound state mass m* fuUfiUs the Strutinsk requirement 
at the extremum 

A™; = . (54) 

However, the pion charge radius calculated according to the non-relativistic formula in p[ is 
about one half of the empirical value. We will confirm this fact, in the renormalized effective 
model while using Eq.(0) and (p.,^ from Eq.(|36|) with A{k) = 1 to compute the pion radius. 

In Figure 2, we present our results for the pion mass as a function of a for rj = 1330 
MeV. Our agreement with the calculation of 0] is within 10% . In Figure 3, we show results 
for our calculation of the ground and excited states masses from Eq. (|53D , for a varying 
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and 7] = 1330 MeV, compared to the renormalized model for the Coulomb plus Dirac-delta 
interaction. In the last case, the bound state masses of the pion ground and excited states 
in s-wave, are found numerically from the zeroes of Eq.(0) with the Green's function of 
the Coulomb-like potential obtained from the solution of the integral scattering equation 
(ppp. In both equations, which define the non-relativistic renormalized model Coulomb plus 
Dirac-Delta interaction, is used A{k) = 1. We have disregarded the effect of the energy 
transfer in of Eq.(|D, as being a short-range effect parameterized by the value of the pion 
mass, which is input in this calculation. In momentum space Eq. (|34|) is given by 

1 1 



oo p2 



poo 

An / dp 
Jo 



A{p) 



m* ^ — Am? — Ap^ — Am? — Ap^ 



2 /•no ^2 



+ 877^ r dq j r dp ^ 



/4(^) ° 

t^{p,q;ml) 



t^{p,q;mf) 



{ml — Am? — 4p^)(m^ — 4m^ — Aq^ 
and s-wave projected T-matrix in Eq.(R3) is 



m*^ — Am^ — 4p^)(m*^ — Am^ — Aq^) 

= ; (55) 



t^{p,q;M^) = J^^dcos{e) <p\T^{M^)\q> ] (56) 



which is the solution of 



Sir^pq JA{p)A{q) 



+ 8m^ dp' p^^^ {p-pT t^{p',q;M^) 

Stt^Vo Ja{p)A{p') P ^ iP + P'y - Am^ - Ap"^ ' 



the momentum space representation of s-wave projection of Eq. (pO|) . 

The agreement between the renormalized model and the smeared delta model is within 
a fraction of percent, which still improves as the weakly bound limit of the ground state is 
approached, as shown in Figure 3. The calculation of is few percent below ours. The 
inclusion of the relativisic phase-space in the renormalized model of the Coulomb plus Dirac- 
delta interaction Eq.(0) make less bound the excited state for a given ground state mass. 
As the pion mass grows the relativistic phase-space effect tends to vanish. 

17 



In the other study performed, a = 0.6904 was kept unchanged, while varying t] between 
350 and 1350 MeV. The results for the excited bound state mass against the ground-state 
mass are shown in Figure 4. The renormalized Coulomb plus Dirac-delta calculation agrees 
within a fraction of percent with the smeared Dirac-delta calculation. The relativistic phase- 
space in the renormalized interaction makes the excited state less bound for a given ground 
state mass. 



C. Effective Meson Model 

The bound state masses of the meson ground and excited states in s-wave, are found nu- 
merically from the zeroes of Eq. (|3^) with the Green's function of the Coulomb-like potential 



obtained from the solution of the integral scattering equation (|20D . The the energy transfer 
in in Eq.(||) was neglected. 



The renormalized strength of the singular interaction, from Eq.(p2D 

=< x|Go(m^)T^(m^)GoK)|x > ; (58) 



with the condition Eq.fpTf) at the physical pion mass is shown in Figure 5 for m = 406 MeV 
and compared to the perturbative calculation 

\n\ml)\,,,, =< x\Goiml)VGoiml)\x > ; (59) 

with V from Eq.(P) and solution of Eq.(^). For values of a below 0.2 the agreement 
between the perturbative and nonperturbative calculation of the renormalized strength is 
quite good, giving confidence to our numerical calculations. 

The wave-function ip„ of the pion from Eq. (|36| ) in momentum space is written as 

m2 -4m2 - V yjA{p) ^ Jo ^ ^^JJ) ^1 - Am? - Ap^ j 



The first-order perturbative pion wave-function is 

1 if 8m r°° dq q (pTgF 

V^7r(p)|pert I A{p) "~ 4m^ — Ap^ I ~'~ 377 JO A{q) p m^ — Am^ — Ap^ ^ ' ^ ^ 



and for a = the pion wave-function has the asymptotic form 



In the pion wave-function eigenfunction of the effective square mass operator of the 
Coulomb plus Yukawa model, Eq.(^), was approximated by the following analytical form: 

¥>:iP) = ,,2 , (63) 

where Af is an arbitrary normalization and the fit is performed for a = 0.6904, t] = 1330 
MeV and m=406 MeV and Pa=515 MeV is adjusted to the numerical solution of Eq . (|53|) . 

The various non-relativistic models of the pion wave-function, with A{k) = l, are plotted 
in Figure 6. In this figure, the comparision between (p.,,, ipn\pert , V'" and ipl are performed. 
The model parameters are a = 0.6904 and m =406 MeV and the normalization is arbitrary. 
The perturbative calculation reproduce ipn{p) for p below m. The asymptotic wave- function 
overestimate (pTr{p), as it should be for a repulsive Coulomb-type interaction, and 
for small p it approaches the nonperturbative eigenfunction. The analytical approximation 
works quite good for momentum up to about 1.5m. More results on the pion wave- function 
are shown in Figure 7, where the results for A{k) = 1 with the full calculation are compared. 
The effect of A{p) diminuishes the magnitude of f{p), as one could antecipate from Eq.(|60|). 
Also we plot ^Pn{p) ioT a = 0.18 and m = 386MeV, which is above the curve of previous 
case due to the decrease of the Coulomb attraction (this parameters fits the empirical pion 
radius). We point also that, for very high momentum, again the asymptotic form dominates, 
as the second term of Eq.(|60|) tends to zero faster than the first one. 

We are going to compute the pion radius in the effective renormalized model of the pion 
using Eq. (|40|) , Eq.(^) and Eq. (^Tf) . In the last two cases, the pion wave-function from ( |60D 
with A{k) = 1 is used. The pion wave-function is known from Eqs. (pSf) and ( |60D from which 
the charge radius including the effect of quark spin is obtained from integration by Gaussian 
quadrature of the difference 

AF^(g2) = F^(g2) _ ^^(0) , (64) 
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obtained from Eq.(^), and the pion radius is calculated from 



(65) 



In the case the quark spin is neglected, two possibilities of calculation of the pion radius are 
considered. One through Eq.(^5D in which only the spin factors are simplified in the limit 
of the quark mass being infinity, while the normalization is defined as in Eg . (|40|) , known 
from the empirical value of and number of colors. The second possibility is to declare 
normalized to one the qq Fock- component of the pion wave-function and use the formula for 
the form-factor in which the quark has no spin, Eq.f^T]). 

The results for the pion charge radius as a function of the strength a of the Coulomb- 
like interaction, Eq. (H), are shown in Figure 8. For = 0.67fm we found a;=0.18 using 
Eq. (^OD . Our calculations obtains the pion charge radius from the relativistic expression (^OD 



which is known to give about twice the non- relativistic radius . The attractive Coulomb 
like interaction increases the radius over the soft-pion limit with a = which is below its 
experimental value. In that sense a consistence is found with the effective theory which 
has an attractive Coulomb-like interaction. The repulsion would be completely inconsistent 
with the pion radius. 

The pion charge radius obtained from first order perturbative calculation of the pion 
wave-function 



^7r(a;, kA_)\ 



pert 



Aik) 

-V-K{^)\pert , (66) 



^ x(l — x) 

is also shown in Figure 8. In consistency with the strength calculation presented in Figure 
5, we observe that for a below 0.2 the perturbative calculations match the nonperturbative 
results. 

The results for the scalar quarks form-factors with from Eq.(pO|) calculated in the 
effective renormalized theory for A{k) = 1 are also shown in Figure 8. The calculation of 
r^r with Eq.(P5D show values above the ones calculated with Eq. (^7D , which are two small 
compared to the empirical value, indicating the importance of the physical normalization 
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of the qq Fock-component of the wave-function using /^r. We also performed the first order 
perturbative calculation of r^^ for the wave-function normalized to 1. The agreement between 
the perturbative nonperturbative calculation is reasonable for a below 0.2 . 

The plot of the pion charge radius against the mass of the first excited-state in shown 
in Figure 9. For decreasing values of a, the pion charge radius diminuishes and consistently 
the excited-state mass increases, i.e., this state becomes less bound. It is clear from this 
figure that to simultaneously fit the radius and the mass of the rho-meson (768 MeV) we 
are obliged to use a different quark mass from the value of 406 MeV. For comparision we 
also show the results from Eg . (|45|) . 

The experimental pion radius of 0.67±0.02 fm is fitted with m„ = m-u = 386 MeV and 
a = 0.18 resulting r^^ = 0.67 fm and the mass of the singlet-2s state of 768 MeV, we remind 
that the pion mass is input in the renormalized model calculation. The singlet-2s excited 
state mass of the uu system is identified with the p meson mass in the present effective 
QCD model. The quark mass is varied to form mesons with one up antiquark together with 
the strange, charm or bottom quarks. The masses of the constituent quarks were within 
the range of 500 to 5000 MeV. The results are shown in Table II and compared to the 
experimental data. The singlet-2s states are identified with the lowest mass vector mesons 
states, since the Dirac-delta interaction is the effective hyperfine interaction, the reason for 
the spliting between the pseudo-scalar and vector mesons. Although, in the singlet channel 
the hyperfine interaction is attractive, which is not valid for the spin one mesons, we believe 
that the Dirac-delta interaction mock up short-range physics beyond that, by taking care of 
the empirical value of the pion mass. 

It is clear that the split between the heavy meson masses will be not adjusted in the 
present calculation for a small. This is reasonable, since these mesons are weakly bound and 
the wave-function spreads out in the region where the confinement potential, not present 
in our model, should be important. In Figure 10, we study the difference of the excited 
and ground state masses as a function of the ground state mass. As we have seen in Table 
II, the difference is underestimated for a = 0.18 above the kaon mass. For a = 0.5 the 
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experimental data is reproduced. This is in fact reasonable if we think that a should on 
average increase with the size scale, indicating the confininig behaviour, which is stronger 
for the heavy mesons since they are less bound than the pion and the kaon. 

VI. CONCLUSION 

The essential development made in this work is the renormalization of the effective QCD- 
inspired Hamiltonian theory with a singularity at zero range and its consequent application 
to the pion and other mesons. The method is an example of the Hamiltonian renormalization 
procedure and it is equivalent to a subtracted equation for the transition matrix. The 
physical renormalization condition is given to the two-particle model at the subtraction 
point and in the cases discussed here it is the ground state binding energy or mass. The 
treatment is shown to be renormalization group invariant, i.e., independent on the arbitrary 
subtraction point. This independence is expressed by a fixed-point Hamiltonian that brings 
the physical input to the theory, the pion mass or ground state binding energy, as well as 
the necessary counterterms that render all the momentum integrations finite. 

First, we have studied in an example, the renormalization method applied to a two- 
body model with a Coulomb plus Dirac-delta interaction, where we have calculated the 
excited state energy for a given ground state binding energy. The results are compared to 
calculations with a Coulomb plus repulsive Yukawa interaction and the renormalized model 
reproduce with less then 10% accuracy the model results of Eq.(^). The success of this 
result drive us to the solution of the renormalization problem of the t|-model 0]. 

The effective mass operator equation in the ||-model Eq.(|l]), has as eigenstate the low- 
est Fock-state component (gg) of the light-cone wave-function of a meson bound system of 
constituent quarks or dressed quarks (not to be confused with the bare quark). This model 
picks out one particular aspect of the gluon exchange between quarks, namely the strong 
attraction of the spin-spin interaction in the singlet channel. Previously [Q], the renormal- 
ization has been carried by first regularizing the Dirac-delta interaction through a Yukawa 
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form and then its parameter was found by the Strutinsky requirement that the mass of the 
excited state is stationary in respect to variation of the regularization parameter. The step 
forward in this work was the use of the renormahzation group invariant approach, in which 
the regularization parameter is not necessary to solve the model. All the short-range physics 
is parametrized by one parameter: the renormalized strength of the Dirac-delta interaction, 
which is determined by the mass of the pion. We showed that, the results for the mass of 
the excited state obtained with the renormalized model and the smeared delta regularized 
model are in agreement within fraction of percent, for the same ground state mass, which 
was varied either by changing the Coulomb interaction intensity or the Yukawa range. The 
concordance still improves as the weakly bound limit of the ground state is approached. 
The effect of the relativisic phase-space in the renormalized model of the Coulomb plus 
Dirac-delta interaction Eg. (p4D makes the excited state, for a given ground state mass, less 
bound. This difference tends to vanish as the pion becomes weakly bound. 

The various models of the qq Fock- component of the pion light-cone wave-function, con- 
sidering A{k)=l and the relativistic phase-space {A{k) from Eq.(^)), were calculated and 
compared with the analytical form and perturbative result. We obtained the reduced wave- 
function ifj, ( pup solution of the renormalized ||-model, (/?7r|pert (0)' analytical form 
(|63|) and the asymptotic form (^). We have used model parameters a = 0.6904 and 
m =406 MeV 0. Although for for p below m all the calculations are reasonably consistent, 
the high momentum tail is dominated by the asymptotic wave- function <f^{p), which over- 
estimate v?7r(p); as it should be for a repulsive Coulomb-type interaction, and for small p it 
approaches the nonperturbative eigenfunction, as well as for very high momentum. The an- 
alytical approximation works quite good for momentum up to about 1.5m, however it does 
not have the asymptotic tail for high momentum. The relativistic phase-space diminuishes 
the magnitude of (p{p) as clearly seen in Eq. (|60|) . 

The calculation of the pion charge radius was performed in the renormalized effective 
Q CD- inspired Hamiltonian theory with Eqs. (BBI) and (EO) and although we have simplified 



the spin dependence in the dynamical equation it is important in the evaluation of the 
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radius. For this purpose we have used an effective pseudo-scalar Lagrangian to construct 
the spin part of the pion wave-function, it gives the absolute normalization of the lowest 
Fock- component of the light-cone wave-function in terms of the weak decay constant /tt, 
the constituent quark mass and number of colors. Turning off the Coulomb-like interaction, 
we retrieved the well knwon result valid in the soft pion limit |T^. We also compared 



the values for r,r obtained with the expressions where the quark spin is neglected and the 
wave-function normalized to 1. We pointed out the necessity of the correct normalization 
of the wave-function of the fj-model according to the /,r value to be able to fit r^, which is 
reminiscent of the fact that /^r and are closely related in the light-front phenomenology 
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The experimental pion radius of 0.67±0.02 fm was fitted with m„ = = 386 MeV and 
a = 0.18 resulting r^^ = 0.67 fm and the mass of the singlet-2s state of 768 MeV, we remind 
that the pion mass is input in the renormalized model calculation. The singlet-2s excited 
state mass of the uu system was identified with the p meson mass in the present effective 
QCD model. We stress that in the singlet channel the hyperfine interaction is attractive, 
which is not valid for the spin one channel of the vector mesons, however, we believe that 
the Dirac-delta interaction mock up short-range physics beyond that, which is brought to 
the model by the empirical value of the pion mass. In essence, without being too naive, 
in order to fit the pion charge radius with sucess and with reasonable parameters, it was 
essential that i) the Coulomb-like interaction be attractive and ii) the normalization of the 
form-factor be consistent with the empirical value. 

In respect to the masses of the ground state of the pseudo-scalar and vector mesons with 
one up antiquark together with the strange, charm or bottom quarks, which were calculated 
as a function of the constituent quark mass, the results were in qualitative agreement with 
the data for a = 0.18. We tried a better fit of this data, since it was clear that the 
split between the heavy meson masses would not be adjusted in the present calculation 
for a small. This in fact, seems reasonable, since these mesons are weakly bound and the 
quarks can be found in a region where the confinement potential, not present in our model, 
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is important. The difference between the vector and pseudo-scalar mesons masses for the 
same qq pair is underestimated for a = 0.18 above the kaon mass. We found that for a = 0.5 
the experimental data is reproduced. This gives us some hope that it is possible to refine 
the tJ,-model to include more physics than initially thought. In regard to this extension, 
it is reasonable to think that a should on average increase with the size scale, indicating 
the confininig behaviour, which we found stronger for the heavy mesons since they are less 
bound than the pion and the kaon. In short, a reasonable description of the physics of the 
pion and other scalar and vector mesons were found, taking into account the simplicity of 
the renormalized effective Light-Cone QCD-inspired theory. 
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APPENDIX A: DERIVATION OF THE RENORMALIZED T-MATRIX 



In this Appendix, the solution of Eq.(0) to find the the renormahzed T-matrix, Eg. (p!9|) , 
is performed in detaiL We want to solve Eq.(p!7D 

Tn{M^) = V + Vl+(y + Vl) G'C\M^)Tn{M^)) • (Al) 

The regular part of the potential, V , is defined by Eq.(|^) and the renormalized singular 
interaction is given by Eg . ([13|) , rewritten as a matrix equation 

vl = \x> Mf^') <x\~\x> Mf^') < xlGUf^'M ; (A2) 

which has the solution 

Vl=\x> vAf^') < X\ , (A3) 

with the function 

(A4) 

The function v-jilfi"^) contains the divergences in the momentum integrals which exactly 
cancels such infinities in Eq. (|Al|) . It is enough for the formal manipulations that will come. 
However, one could equally well introduces a cutoff in Eqs. (|A1|) and ([A^), and performs 
the limit of the cutoff going to infinity just after the solution of Eg .([All), in which all the 
necessary cancellations happens and the limit is finite. 
Next, Eq. (|Al|) is rewritten as: 

(l - VGi-'\M')) Tn{M') = V + Vl + VlG'^^\M^)Tn{M^)) , (A5) 

and inverting the operator in the left-hand using the regular T-matrix, T^(M^), solution of 
Eq.(pOl), one has 

Tn{M'^) = T^iM^) 

+ (l + T^(M2)G^+)(M2)) \x > vnifi') < x\ {G^^\M^)Tn{M^) + l) . (A6) 
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The "bra" function, < x\Go'\M'^)Tn{M'^) ^ has to be calculated in order to find the 
renormalized T-matrix. We multiply Eg . ([A^) by < x\G-^o^\m'^) on both sides, and solving 
it we get 

<X\G, {M)Tn{M)- l-<x\GyM{M^)\x>vn{^^') ' ^^^^ 

where regular potential resolvent is 

G^^+\M^) = G^^\m^) + G^^\m^)T^ {M^)G^^\m^) . (A8) 

The "bra" function of Eq. (|A7|) is introduced back in Eg . (|A6|) , and with a little algebra 
one finds 

(l + T^{M^)Gi^\M^)) \x >< x\ (g[,+^(M2)T^(M2) + 1 
v^\^I^)-<x\G''MiM^)\x> 



which, after introducing Eg. ( |A4D and the explicit form of the resolvent of Eg. (|A8| ) , results 
in the renormalized T-matrix of Eg. (|19D . 
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TABLES 



r/ e(i) [C-Y] e(2) [^.y] e^^) [C-6] 

0.1 0.1109 0.06781 0.06237 

1 0.5119 0.1813 0.1736 

10 0.9495 0.2449 0.2439 

TABLE L Results for S-wave binding energies of the non-relativistic model of Eq.(^), our 

calculation with N = 100 [C-Y], compared to the effective Coulomb plus Dirac-delta interaction 
[C-6]. 



qq 


mi 


m.2 




b,exp 




b,exp 


ud 


386 


386 


140 


TT^ : 140 


768 


p : 768 


us 


386 


500 


511 


: 494 


882 


K* : 892 


uc 


386 


1500 


1852 


L>° : 1865 


1882 


L>*° : 2007 


ub 


386 


5000 


5375 


: 5279 


5383 


B* : 5325 



TABLE II. Results for S-wave Meson Masses: singlet-ls (M^^"^) and singlet-2s iM^f^) and 
experimental values {Mj^^^^^^^). All masses in MeV. a =0.18. 
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Fig. 1 Excited state binding energy (e^^-*) as a function of the ground state binding energy (e^^-*) 
for the non-relativistic model. Attractive Coulomb plus repulsive Yukawa non-relativistic 
model of the Ref. ||l^ (full square). Our numerical calculation of Eq.(|48|) is given by the 
solid line, and the solution of the effective model Eq.(^) (Coulomb plus Dirac-Delta) is given 
by the dashed line. 
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Fig. 2 Pion mass as a function of the strength of the Coulomb interaction a. Non-relativistic 
model for Coulomb plus Yukawa interaction, m = 406 MeV, solid line. Yukawa range 
parameter 77= 1330 MeV The full square is the calculation of Ref. Q. 
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Fig. 3 Mass of the excited state as a function of the ground state pion mass. Non-relativistic 
model for Coulomb plus Yukawa interaction, m = 406 MeV, solid line. Strength of the 
Coulomb potential within < a < 0.7 and Yukawa range parameter r/=1330 MeV [||. The 
full square is the calculation of Ref. |^. The non-relativistic effective model (Coulomb plus 
Dirac-Delta) is given by the full dots. The relativistic effective model for the Coulomb plus 
Dirac-Delta interaction is given by the dashed-line. 
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Fig. 4 Mass of the excited state as a function of the ground state pion mass. Non-relativistic 
model for Coulomb plus Yukawa interaction, m = 406 MeV, solid line. Strength of the 
Coulomb potential a = 0.6904 and Yukawa range parameter within < rj < 1350 MeV 0. 
The full square is the calculation of Ref. j^. The non-relativistic effective model (Coulomb 
plus Dirac-Delta) is given by the full dots. The relativistic effective model for the Coulomb 
plus Dirac-Delta interaction is given by the dashed line. 
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Fig. 5 Inverse renormalized strength X^{m^) in units of as a function of the Coulomb intensity 

potential parameter a for a pion mass of 140 MeV. Nonperturbative calculation (solid line) 
and first-order perturbative calculation (dashed line), see text for the explanation. 
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Fig. 6 Pion wave-function ((/?) witii arbitrary normalization as a function of momentum in units 
of quark mass. Calculations performed with A{k) = 1. Pion model wave-function model of 
(|60| ) for a = 0.6904 and m = 406 MeV (solid line), first-order perturbative calculation from 
Eq.(^) (short-dashed line), asymptotic form Eq. ( |62|) (dotted line) and fit from Eq.(|63|) 
(dashed line). 
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Fig. 7 Pion wave-function (ip) as a function of momentum in units of quark mass. Pion model for 
a = 0.6904 and m = 406 MeV, calculation with Eq.(|60|) (dashed line); considering A{k) = 1 
(solid line). Results for a = 0.18 and m = 386 MeV (dotted line). 
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Fig. 8 Pion charge radius as a function of the strength a of the Coulomb potential. Quark mass 
of 406 MeV. Results for Eq.(]40|): nonperturbative calculation of the wave function (solid 
line) and first order perturbative calculation (dashed line). Results for Eq.(|45|) and wave- 
function obtained with A{W) = 1 (solid line with triangles). Results for Eq. (|47| ) with the 
wave function obtained with ^(fe) = 1: nonperturbative calculation of the wave function 
(solid line with dots) and first order perturbative calculation (dashed line with dots). 
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Fig. 9 Pion charge radius as a function of the mass of the excited qq state for a constrained pion 
mass of 140 MeV. The strength a is within the range of Fig. 6 and the quark mass is 406 
MeV. Results for Eq.(||) (soHd Hue) and for Eq.(|5l) with wave-function obtained considering 
A{k) = 1 (soUd Hue with triangles). 
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Fig. 10 The difference between tlie masses of tlie excited and ground state of the qq system as 
a funtion of tlic ground state mass in the effective relativistic model. The quark mass is 
varied in the range between 386 and 5000 MeV. The other quark mass is fixed to 386 MeV. 
Calculations with a = 0.18 (solid line) and a = 0.5 (dashed line). The experimental values 
from Table II are given by the full circles. 



39 



